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Abstract 

The magnetization and the de Haas-van Alphen oscillations of 
Bloch electrons are calculated near commensurate magnetic fluxes. 
Two phases that appear in the quantization of mixed systems — the 
Berry's phase and a phase first discovered by Wilkinson — play a key 
role in the theory. 

The magnetization of a free electron gas was calculated by L. Landau 
in 1930 in the early days of quantum mechanics [1]. Considerable efforts 
have since been devoted to extending Landau results to Bloch electrons, i.e., 
in the presence of periodic background potential. Most of the efforts and 
progress made was in the region of weak magnetic fields [2, 3] where the flux 
$ through a unit cell is small. This is adequate for most solid state applica- 
tions. There is, however, also considerable interest in a better understanding 
of phenomena that have to do with commensuration in condensed matter 
physics [4, 5, 6, 7, 8]. This is the case when $ is close to a rational num- 
ber. The magnetization of Bloch electrons near rational fluxes, other than 
$ = 0, remained an open challenge which we solve here. The difficulty lies 
in the delicate spectral properties resulting from commensuration [5]. The 
Hofstadter model is a basic model for a system where commensuration plays 
a role. It is also a basic paradigm for a quantum system with fractal spectra 
[6], anomalous quantum transport [7] and the (integer) quantum Hall effect 
[9]. 

The problem of magnetization near fractional flux becomes tractable by 
an idea that goes back to M. Wilkinson [8]. Namely, that near a rational 
flux the Hamiltonian can be understood as the semiclassical quantization of a 
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mixed system: In mixed systems some, but not all, degrees of freedom may be 
treated semiclassically. As a consequence the "classical Hamiltonian" is ma- 
trix (or operator) valued. Pauli and Dirac equations for a spinning electron 
in a slowly varying potential, the Born-Oppenheimer theory of molecules, 
and the Hofstadter model near rational flux [8] are examples of mixed sys- 
tems. In the Hofstadter model the role of the Planck constant is played by 
the deviation from a nearby rational 

h = $-- . (1) 

Littlejohn and Flynn [10] developed an elegant geometric formalism for 
the quantization of mixed systems. They show that in order h the quantiza- 
tion of mixed systems gives rise to two phases: One is the Berry's phase [1 1] 
and the other is a phase that is sometimes known as the "no-name phase" [12] 
and sometimes as the Wilkinson- Rammal (WR) phase [13]. For both phases 
to appear, the "classical Hamiltonian" must have non-trivial commutation 
properties in both coordinates and momenta. These phases play a central 
role in determining the magnetization, see e.g. Eq. (25) below. 

Our results are closely related to recent progress made in the semiclassi- 
cal dynamics of Bloch electrons under slowly varying electric and magnetic 
fields [13]. When one goes beyond the leading order expressed by Peierls 
substitution [2, 13] one finds that the Berry's phase and the WR phase play 
a role in the dynamics. This lead [13] to identify the WR term with the mag- 
netization of a wave packet. Although related, the notions of wave packet 
vs. thermodynamic magnetization, expressed in Eq. (25), are distinct; for 
example, wave packet magnetization is not defined in the gaps, while the 
thermodynamic magnetization and, of course, the de Haas van Alphen oscil- 
lations, have a non-trivial dependence on the chemical potential also in the 
gaps. 

Fig. 1 shows the zero temperature magnetization at $ = tit-. The 

complexity of the magnetization is due to the multiplicity of scales: The big 
scale is determined by the denominator q = 3, and the small scale by h. On 
the small scale one sees the rapid de Haas-van Alphen oscillations. On the 
big scales one sees continuous features: the linear pieces in the (big) gaps, 
the envelopes of the de Haas-van Alphen oscillations and their mean. Our 
theory of magnetization accounts for all these features. 

We show that while the amplitude of the oscillation is determined solely 
by the leading terms in the semiclassical expansion, the mean magnetization 
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requires knowledge of the terms beyond leading order, thus depending on 
the fine details of the spectrum. Intriguingly, it is the latter quantity which 
is stable against perturbations. Finite temperatures larger than the typical 
eigenvalue spacing wash out the de Haas- van Alphen oscillations of the small 
scale but leave intact the mean magnetization. Semiclassical approximations 
that retain only the leading order yield no magnetization at all at finite 
temperatures. 

Let us start by describing the semiclassical quantization of mixed sys- 
tems [10]. The "classical Hamiltonian" , Ti(x, k), is a Hermitian matrix which 
depends on x and k. We shall denote by £j(x, k) its j-th band of eigenvalues 
and by \uj{x, k)) the corresponding eigenvectors. We shall also assume that 
bands do not cross 1 . The corresponding quantum Hamiltonian is Ti(x, k) 
with [x, k] — ih. 

Let S(E; H, p) denote the classical action associated with a closed orbit 
of energy E of a classical Hamiltonian function H, with phase space area 
form p. Since phase space is two dimensional the action is the area enclosed 
by the orbit. Note that H, unlike 7i, is a scalar valued function. The Bohr- 
Sommerfeld quantization rule in mixed systems says that the semiclassical 
approximation of the eigenvalue E n is given by 

S(E n ;H h ,p h ) = h{n + lm ) , neZ (2) 

where 7 m , is the Maslov index of the orbit [14] . has an expansion in powers 
of h, H h = if(°) + HH^ + ■■■. Peierls substitution sets Hf\x, k) = Ej(x, k), 
and the next order is [10] 

Hj 1 \x, k) = Im (d x uj\ (H — Sj) \ d k Uj) (x, k) . (3) 

The expansion pn, = p^ + hp^ + • • • begins with the canonical form p^ = 
dk A dx and the subleading term is the Berry curvature form 

p^P = 2uj(x, k) dk A dx , u(x, k) = Im (d x Uj\ dkUj)(x, k) . (4) 

This formulation is manifestly gauge invariant (independent of the choice 
of phases for \uj)) and preserves the symmetry properties of 7i, which is 
useful when one wants to correctly count the dimension of the Hilbert space 
of the quantized operator, as we now proceed to explain. 

1 In the Hofstadtcr model, this is guaranteed by Chambers relation. 
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Suppose that Ti is periodic in both x and k up to gauge transformations, 
and hence describes (classical) motion on a phase space torus T. Hj°\ 

and rfp are all well-defined functions on T. The Chern number of the j-th 
band is the integer 

(5) 



7T Jr 



ujj dk A dx 



It follows from Eqs. (2,4) that the dimension of the Hilbert space associated 
with the j-th band is 

X + C7 '"' (6) 

where |T| denotes the area of T. Since the dimension of the Hilbert space is 
necessarily an integer, quantization on the torus is possible only for certain 
values of h. Eq. (6) goes beyond the classical Weyl law which only determines 
the leading, \T\/h, behavior. The Chern numbers shift states between the 
spaces of different bands since J^Cj = [15] . 

The H corrections to the action can be moved from the left hand side of 
the Bohr-Sommerfeld relation to the right hand side, where they acquire an 
interpretation as two additional phase: 

S J (E)=h(n + lm + T 3 (E)), -2rrT J (E)= 1 b J (E)+ 1 J r (E), (7) 

where Sj(E) = S(E; £j, dk A dx). rf(E) is the Berry's phase [11], 



(8) 



7 j(£) = 2 J ojj 6(E - ej) dk A dx, 
and 7j or (£') is the Wilkinson- Rammal phase 



~fJ r (E) = - J Hf ] 5{E- Ej) dk A dx. 



(9) 



It is noteworthy that the WR phase need not vanish at band edges. 

Let us now recall some basic facts about the Hofstadter model [5, 8]. 
When the magnetic flux $ is a rational number p/q, the model is represented 
by [8] 

H(x, k) = e 2nix U + e 2nik V + h.c. (10) 
where U and V are the q x q matrices 
/0 1 ••• 0\ /l 

/ .'' , v = 



\1 



1 

0/ 



e 2-Kip/q 



\ 



e 2wi(q-l)p/q j 
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The magnetic bands of the Hofstadter Hamiltonian at $ = p/q are given by 
Sj(x, k) on the Brillouin zone 

BZ = {x,k\0 < x < 1, < k < 1/q}. (12) 

Evidently, H(x, k) is periodic with period 1 in both variables. Moreover, 
H is periodic with smaller periods up to unitary transformations: 

H(x, k) = TH(x + 1/q, k)T^ = GH(x, k + l/q)G ] (13) 

G is a gauge transformation (a diagonal unitary) and T a shift. This makes 
the band dispersion functions Sj (x, k) periodic with periods 1/q in each vari- 
able and with q periods in T = BZ. 

The spectrum of the Hofstadter model for other values of <£> is obtained 
by setting [x, k] — %h in H. with h given by Eq. (1). BZ is the minimal 
torus on which 7i may be quantized. The dimension of the Hilbert space 
associated with it is then — Cj. Since the band function are g-periodic on 
BZ, the number of distinct eigenvalues is of order l/(q 2 h). The semiclassical 
approximation is valid provided this number is large i.e. q 2 h <C 1. 

We now turn to the magnetization of the model. Recall that the Hof- 
stadter model approximates the Schrodinger equation in two dual limits: 
When the magnetic field is weak relative to the periodic potential and also 
in the opposite limit where the magnetic field dominates all other interac- 
tions. The two limits have related but different thermodynamics. For the 
sake of concreteness we shall consider the tight-binding interpretation. The 
magnetization of the "split Landau level" follows from the duality transfor- 
mation of [16]. 

The thermodynamic potential per lattice site of the Hofstadter model for 
rational flux and zero temperature is [16] 

n(fi,$) = J dkAdxTr^-T^j (14) 

where x + = x9(x). When jj, is in a spectral gap, the thermodynamic potential 
can be written as a sum of the potential of the occupied bands, fl = J2 < f2<. 
The magnetization per unit area m is 2 



2 To translate the magnetization to ordinary units one needs to divide our dimensionless 
magnetization by the unit of quantum flux. 
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Figure 1: The de Haas van-Alp hen oscialltions of the magnetization as a 
function of the chemical potential in the Hofstadter model for flux $ = 
40/121, compared with the limiting envelopes and their averages, and the 
gap magnetization for $ — > 1/3. The two insets are details of the left and 
center envelopes. 

The magnetization in the gaps can likewise be expressed as a sum of the 
magnetization of the occupied bands 



where the magnetization of a full band is, as we shall see below, (see also [16]): 



The term proportional to /i is the Chern number of the band. It follows that 
the magnetization as a function of // has quantized slopes in the gaps. 

The envelope of the de Haas- van Alphen oscillations of j-th band, as we 
shall show below, is given by 




(16) 



< 




(17) 



Lj((jt, $) = (m)j(fi, $) ± SrrijifjL, $). 



(18) 
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(m)j(/x, $), is the natural restriction of Eq. (17) to a partially filled band, 
i.e., 

(m)j{fji, $) = -— / 9(n - 8^(2(11 - ej)uj - HfA dk A dx. (19) 

27T J BZ V / 

It describes the mean value of the magnetization, averaged over the de Haas- 
van Alphen oscillations. 

The width of the envelope is given in terms of the classical action associ- 
ated with ef 

**o.) = ffjg- (20) 

^ S'(fx) is proportional to the density of states. Since the density of states 
in two dimensions diverges logarithmically near the separatrix, the width 
5m shrinks to zero logarithmically there. Near the bottom of the band Sj 
vanishes linearly while the density of states approaches a positive value. This 
shows that 5m vanishes linearly at band edges. These properties characterize 
the universal lip-like shape of the envelopes. 

We conclude with an outline of the derivation of Eqs. (18,19,20). Consider 
the zero temperature thermodynamic potential associated with one fixed 
band j. It follows from the Chambers relation [17] and the square symmetry 
of the Hofstadter Hamiltonian that for all energies except the separatrix, the 
level sets of Ej(x, k) are deformed circles, and therefore 7 m = 1/2. All spectral 
quantities below refer to the same band, and we may therefore suppress the 
index j without risk of confusion. 

Suppose that /i is such that n spectral points of the split j-th band are 
occupied. Recall that each spectral point is g-fold degenerate and that, by 
Eq. (6), nq < ^ + C. Suppose for definiteness that /i is below the seperatrix. 
By the Bohr-Sommerfeld rule the thermodynamic potential is (to order h 2 ) 

n($,fi)=qh - S ^ (( £ + V2 + T(E e ))h^j , (21) 

The overall factor qh comes from the degeneracy per unit area of each eigen- 
value. Approximating the sum with the second Euler-Maclaurin sum formula 
gives (again to order h 2 ) 

pnh 

n($,fj.) =qfj,nh-q dx S~ l (x + hriS'^x))) . (22) 
Jo 
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Taking derivative with respect to $ is the same as taking derivative with 
respect to h. The magnetization is therefore given (to order K) by 

m(Q,(i) = -qn(n- ^((n + 1») +q J\(E)dE 

= qn(^ En+ 2 En ^ -f^j+g fr(E)dE, (23) 

where the E'-integration is over the j-th band energies below \x. 

The first term describes the de Haas-van Alphen oscillations: It vanishes 
in the middle of each spectral gap [E n -i,E n ], and reaches its maximum 
magnitude at the band edges. 8m{n) is half the variation of m across a 
spectral gap. Therefore 

M<M=?" 2 ^25V) (24) 

The first term in Eq. (23) has zero average over the gap, while the second 
term is nearly constant. The mean magnetization is therefore 

(m) (<£>,//) = - J- y" ^(E) + r r (E)) dE (25) 

The Berry's phase contributes 

q f 1 \E)dE = 2 f udkAdx 6(E - e)dE 
J Jbz J 

= 2 / (n - e)u9(n - e) dk A dx (26) 

JBZ 

The WR-phase contributes 

q [% wr (E)dE = - f dkAdxH {1) f^dE5(E-e) 
J Jbz J 

= - [ H (1) 6{E - e) dk A dx (27) 

J BZ 

Together, they add up to give Eq. (19). 

Finally, let us present a streamlined derivation of the rules for band split- 
ting [8]. Consider for example, 

$ = r = T = - + K h = , (28) 

q ~ $ qn-1 q q(qn - 1) 



with q odd and n even. We demonstrate the following splitting rule: Of 
the q bands associated with the flux 1/q, the center band splits into n — 1 
subbands and the rest into n subbands, together accounting for the qn — 
1 band associated with the flux $. Clearly, this should follow from the 
dimension formula Eq. (6), which requires the additional input of the Chern 
numbers. The Diophantine equation of [9] at flux 1/q bands implies that the 
Chern number of the center band is 1 — q, and all other bands have Chern 
number 1. Recalling that the area of the Brillouin zone is 1/q, Eq. (6) gives 
that the center band splits into q{n — 1) levels and the other bands split 
into qn levels each. We recall also that the band dispersion functions have 
q periods in the Brillouin zone, and therefore there are only qn — 1 distinct 
levels, each g-fold degenerate 3 . This example illustrates the algorithm which 
generates the hierarchical structure of the Hofstadter butterfly. 
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